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Abstract 

We review an exact analytical resolution method for general one- 
dimensional (ID) quantal anharmonic oscillators: stationary Schrodinger 
equations with polynomial potentials. It is an exact form of WKB 
treatment involving spectral (usual) vs "classical" (newer) zeta-regula- 
risations in parallel. The central results are a set of Bohr-Sommerfeld- 
like but exact quantisation conditions, directly drawn from Wronskian 
identities, and appearing to extend Bethe-Ansatz formulae of inte- 
grable systems. Such exact quantisation conditions do not just select 
the eigenvalues; some evaluate the spectral determinants, and others 
the wavefunctions, for the spectral parameter in general position. 

It is a great honour and pleasure to write in this volume for Professor 
Dowker as a tribute to his prominent contributions to mathematical physics. 
Because of our common interests toward functional (i.e., zeta- regularised) 
determinants in spectral theory, we stay in this area. However, instead of 
the more usual applications in geometry (e.g., [18, 9]), quantum fields (e.g., 
[12, 14]) and analytic number theory (e.g., [10, 26]), here we propose a foray 
into analysis, with an exact treatment for Sturm-Liouville equations — or 
ID stationary Schrodinger equations in quantum mechanics. 

We specifically review our exact- WKB resolution of that problem, trying 
to draw the main lines in a global perspective, which is blurred by many 
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technicalities in our earlier works. Therefore we make no attempt at com- 
pleteness here and may discard even important details for which we refer to 
earlier presentations ([20, 21, 23, 25] on the technical side; [22, 24] on the 
survey side). Likewise, the bibliography here is partial and largely meant to 
refer to those in our earher works, most developed in [21, 24, 25]. Inversely, 
two points previously scattered among several articles are revisited here in 
greater detail: the classical side of zeta-regularisation in Sect. 2, and the 
proof of the basic exact- WKB identities (32) in Sect. 3 and the Appendix. 



1 General setting 

We review an analytical resolution method, through an exact- WKB approach 
assisted by zeta-regularisation, for the stationary ID Schrodinger equation 
(or Sturm-Liouville problem) 

with V{q) (the potential) a polynomial function, here normalised as 

V{q) ^+q^+vi g^-i + • • • + Vn-i q {N>1). (2) 
Key parameters will be 

N (degree), ^ + ;^ (order), (f =^ (symmetry angle), (3) 

plus the "number of conjugates" L to be determined in (5). Often we will 
simplify results by excluding N — 2, which is most singular as special case 
within the general- A?" theory. 

Other quantum-mechanical notations will also help: H (Schrodinger op- 
erator), E (energy variable), condensing (1) to 

Hij = Eij, with H=-$^ + V{q), E = -X, (4) 

dq'^ 

but this continues to mean the bare differential equation (1): is not re- 
stricted a priori to eigenvalues, nor 'i/j to some operator domain. 

Being motivated by self-adjoint bound-state calculations in quantum me- 
chanics, we envision both q and V as real, or -u =^ (vi, . . . , vn-i) £ M-^~^. For 
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such spectral problems the best setting is a real half-line, say q G [0, +00), 
with ultimately a Dirichlet or Neumann boundary condition at the endpoint 
g = 0. Still, our exact resolution will involve 

- the whole complex g-plane; 

- all "conjugate equations" of the form (1) but complex, that restore (1) 
under some analytic dilation q i-> e'^g; [15, §2.7] there are finitely many, 
precisely 



distinct conjugates : (5) 



N + 2 in general 

+ 1 for an even polynomial V{q) 

V^%)^^e-''^V{e-''^/'q), X^'^^^e-^'^X for £ = 0, 1, . . . , L - 1 (mod L), 

(6) 

with the angle <^ = 47r/(A^'-|-2)as announced in (3). 

Traditionally, the ordinary differential equation (1) only admitted asymp- 
totic (hence approximate) solutions in analytic form: WKB expressions, valid 
for g or A large. Then, parts of WKB theory have been turned fully exact 
[6, 3, 15, 17, 1]. What we will outline here is a global exact-WKB treat- 
ment for (1), which just falls short of self-sufficiency: curiously, it still needs 
explicit inputs of prior asymptotic- WKB results; hence we review these first. 



2 Asymptotic ingredients 
2.1 Known spectral properties 

To streamline this part, we phrase it for even polynomials V{q). (The non- 
even case adds minor, but tedious, comphcations. [21]) Then, the Schrodinger 
operator H over all of R, with V{q) f +00 at both ends q — >■ ±00, has a 

discrete countable £'-spectrum S =^ {Ek}k=o,i,... with Ek t +00. By parity- 
symmetry, moreover, '= {-Efcjfe even is the Neumann spectrum for (4) over 
the half-line q e [0, +00), while S~ =^ {Ek}k odd is the Dirichlet spectrum. 



2.1.1 Semiclassical asymptotics 

The spectral semiclassical regime is A; — )■ +00, for which the eigenvalues Ek 
asymptotically obey a Bohr-Sommerfeld condition shaped by the classical 
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~ + I for integer k — >■ +00; (7) 



Hamiltonian + V{q) over the phase space R^, as [5, §10.5] 

/ pdq 

Jp'^+v{q)=E 27r 

the leading power behaviour in £' ^ +00 is 

pdg L 771U -j-i- ^ , ^ /„JL / P*^^ 



2 

E=Ek 



i ^ ~ KE'' with = 1 + i- (and 6^ = / 



(8) 



introducing an essential parameter, the order e (g, |]- However, a stronger 
result needs to be known: that the complete Bohr-Sommerfeld formula, 
which incorporates the semiclassical corrections to (7) to all orders [5, §10.7], 
admits a complete E +00 expansion as a descending power scries, 

^baEk"r^k + ^ for integer /c +00, A ^ {/i, /i- ^, /i- ^, . . .} {9) 

(each ba is a polynomial expression in the {vj}j<(ij,-a)N)- 
2.1.2 Spectral functions 

The form (9) implies the following statements, issued for real A > — min V 
initially (other A are to be reached by analytical continuation). 
- The pair of series 



Z^{s, A) = J2 + A)-' for Res> n (10) 

k 



odd 



converge, defining analytic functions Z"^: (generalised) spectral zeta- functions. 

- The summations in (10) can be handled by Euler-Maclaurin formulae, 
whereby continue to meromorphic functions in the whole complex s-plane, 
pole-free (i.e., regular) at s — 0. [20, §1.2] 

Besides, (term by term) all sums of the form (10) obey the general identity 

dxZ{s,\) = -sZ{s + l,\). (11) 

Next, the regularity of Z"^ at s = allows to specify spectral determi- 
nants, not just formally as 

D±(A) = " Yl i^+Ek) " (severely divergent), (12) 

. even 
" odd 



but rigorously, as well-defined entire functions of A, through 



logD±(A) 



def 



dsZ^{s,X)s=o 



{zeta-regularisation) . 



(13) 



In (13), the outputs of Euler-Maclaurin summation are limit formulae 
which we call structure formulae (akin to Hadamard products, but with no 
free constants left). 



with parity preserved: k, K in (14) are even for , odd for D~ . 

Finally, \ogD^ have (term-by-term computable) A — )■ +oo expansions 
called generalised Stirling expansions [13], which get constrained by the zeta- 
regularisation to a very specific form: 



(A") = A" save for a e N : (A°) = log A, (A^) = A(logA - 1), ... (16) 

(...(A") = A"(logA — Hn), Hn — the harmonic numbers), but only one 
implication will matter here: that no stand-alone powers hn\^ {n G N), 
including additive constants (boX^), can appear. We dub the restricted form 
(15) "canonical"; it actually stems more directly from a (more general) heat- 
trace expansion than from (9) [21, eqs. (8)]. 

At N — 2 (the harmonic-oscillator case), classic F- function results are 
recovered. Either spectrum, say S, essentially amounts to the set N, for which 
/i = 1, the generalised zeta function is Hurwitz's Ci^^^^) = X]fclo(^ ^)~*) 
and the zeta-regularised determinant is e"^*^^^''^^ = \/27r/F(A). [11, vol. I 
eq. 1.10(10)] Then (14), resp. (15), restore the Euler limit formula, resp. the 
Stirling expansion, for — logF(A). [11, vol. I eqs. 1-1(2), resp. 1.18(1)] 





over the set A of (9), with 



(15) 



5 



2.2 Wave asymptotics 

2.2.1 Known WKB results 

The differential equation (f) on the half-line [0, +00) has its solutions if) 
asymptotically described (for q — >■ +00) by WKB formulae, 

m oc n,(g)-V2 exp|± fu^mq^, Mq)'^ (V{q)+Xy/' ; (17) 

i.e., inA(?) is the momentum in the classically forbidden region; primitives 
(as in the exponent) are the action functions; Il\{q) > for q — > +00 fixes 
the branch. Differentiations of (17) (with respect to q. A) arc also valid. 
An explicit q ^ +00 form for stems from the classical expansion 

[V(q) + Xr+y'r.J2^a{s;v,X)q''-''\ a=f, f-1, f-2, ... (18) 

c 

(in which the come out as polynomials in v and s), as 

i){q) oc g-^/4±^-i(^^)exp [±^fi„_^{v) ^| (19) 

^ {a>0} ^ 

[15, Thm 6.1(ii)] where ^„_\{y) ^= /3ct-i(0; A) is independent of A for cr > 
when N ^2. 

2.2.2 "Cleissical zeta-regulEirisation" 

The overall constant factor imphed in (19) depends on the primitive of 
±Il\{q) selected within (17). Against common practice and for later sim- 
plicity, here we opt for an "improper action", formally meant as [21] 

r+00 

Ix{q) = Ilx{q)dq" (a highly divergent integral). (20) 

Jq 

A natural regularisation for (20) would be the analytical continuation of 

r+oo 

h{q, s) / {V{q) + X)'/'-' dq (Re s > /x) (21) 

Jq 
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to s = 0, but the singular structure of Ix{q, s), which follows from the expan- 
sion (18) term by term, shows (for even N) a generic simple pole at s = of 
residue N~^l3-i{v); just for definiteness, (3-i{v) is the (finite) sum [20] 

(22) 

Thus (even though the Kronecker S in (22) can often suppress the whole 
sum), a fully general specification of Ix{q) must accommodate P^i{v) ^ 0; 
this we pursue again at g = and for an even V{q) (still for simplicity only). 

We then argue that (20) (at q = 0) can be fully defined in complete 
classical parallelism with the zeta-regularisation of (12) to logD(X) in the 
quantal setting, when D{X) = det{H + A) over the whole real line or D{X) = 
D^{\)D~{X), all under V{q) even. Indeed, the classical limits of the Schrodinger 
operator H and of the trace operation Tr arc the Hamiltonian function 
+ V{q) and the phase-space integration (27r)^^ J dpdq respectively. Hence 
the classical analogue of the zeta-function Z{s, A) = Tr (^ -|- A)~^ has to be 

Zciis, A) = / dgdp (p2 + Viq) + A)-^ (Re s > //) 

and the last expression yields a regular continuation at s = (where the 
simple pole of Ix{q, s) gets cancelled by r(s)~^). We can thus define 

logL'ei(A) -dsZais, \)s=o. (24) 

Now when Ix is regular at s = 0, — |(9sZci(s, A)s=o evaluated through (23) 
coincides with lx{0,0) (= jQ°°Ilx{q) dq formally). Inversely then, 

r+oc 

/ Ux{q)dq = llogD^i{X) (always finite) (25) 
Jo 

extends the improper action to the singular- 7^ case as well. Equivalently: in 
parallel to the quantal case (although has a branch cut instead of being 
entire) logDci has a canonical large-A expansion; now only one primitive of 
Jq°° dx[lix{Q)] dg can have that property; it then specifies J^^°° lix{Q) dg. 
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Improper actions benefit WKB theory in many other ways, [25] and they 
are no less accessible than conventional action integrals: e.g., 

^ ^ _ r(l + Jf)ji-l - ^) ^1+^ (^_,2) (26) 

Jo 2V7r 

= -|A(logA-l) (TV = 2) (27) 

r+oo 

/ {q^ + vq^ + Xy^^ dq = (case A > 0) 

Jo 

liv + 2y/Xy/^[2VXK{k) - vE{k)], k = > 2VA) 

I AV4[(2x/A + v)K{k) - 2vE{k)], k {v < 2^) 

(28) 

where £"(•), K{-) are the usual complete elliptic integrals. [11, Vol. II, §13.8] 
The calculation as above but for the ratio / D~ in place of the prod- 
uct D^D~ is now regular, and its output, [Z^^/D^](A) = nA(0), completes 
the results we need to rewrite the WKB formulae suggestively for later anal- 
ogy. Together with (25), and with translation covariance restored (thanks 
to Jj-Tlx{q)dq being canonical for finite intervals X), all of that converts 
the canonically normalised recessive WKB-solution (17) plus its g-derivative 
partner, namely 

V'wkb(?) = +n,(?)-V^ exp/+°°nA(g)dg, 
[^VKB(g) = -Iix{qr"' exp/+'-n,(g)dg, 
to the form (also valid for non-even V) 

V.WKB(g) = ^-ci(A), WWB{q) ^ -D+JX), (30) 

where D^^[(A) are the "classical determinants" defined like -D^(A) but over 
the half-line [q,+oo) in place of [0, -|-oo). 



(29) 



3 Exact analysis of the Schrodinger equation 

The differential equation (1) is known to admit exact solutions ip{q) which 
are recessive, i.e., decaying for q — > -|-oo, according to the asymptotic WKB 
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form (17). We specifically select the exact solution ip\{q) that is canonically 
recessive: i.e., asymptotic to our normalised (decaying) WKB form (29), 

r+oo 

i^xiQ) ^ V'WKB(g) = nA(g)"^/^ exp / Iix{q)dq for q +00. (31) 

Jq 

3.1 Basic exact- WKB identities 

The exact solution ipx, given the normalisation (31) for q — >■ +00, is then 
expressed at finite q by the exact identities [21] 

Uq) = ^,-(A), i^'M ^ -D+iX), (32) 

in terms of the zcta-rcgulariscd Dirichlet (— ), resp. Neumann (+), spectral 
determinants over the half-line [q, +00). 

The crucial feature of (32) is the unseen one! No fudge factor whatsoever 
in-between the solution (^a,V'a) the spectral determinants D^, thanks 
to careful normalisations, is what will allow for exact solution algorithms. 

In view of the alternative form (30) exhibited by the usual asymptotic- 
WKB solutions, we see (32) as exact-WKB formulae for the exact solution ipx. 
Clearly these do not appear as explicit as the asymptotic- WKB solutions, 
expressible by quadratures as in (29). Nevertheless, our main point is this: 
(32). too, will end up "solving" the Schrodinger problem analytically, and in 
a more roundabout but now exact, manner. 

Proof of (32) (adapted from [17, Apps. A & D]). It has two steps: 1) 
prove the first logarithmic derivatives of the identities (32) with respect to A; 
2) do one controlled integration thereupon to attain (the logarithm of) (32) 
itself. (All that for N > 2; here we skip the particular cases = 1 and 2, 
which need two differentiations/integrations.) We denote' = ^ and, without 
loss of generality, set g = again. 

1) Applying to (13) the basic functional relation (11) tuned to s — 0, 
[9sZ(s, A)s=o] — —Z{1,X), we find the logarithmic derivative of (32) with 
respect to A, at g = 0, to be 

V'a(0)/Va(0) = Z-{1, A), V'UO)M(O) = Z+{1, A). (33) 

We then refer to the Appendix for a trace-formula proof of (33). 

2) One integration upon (33) with respect to A gives 

logV'A(O) = logD-{X) + C-, log[-V'l(0)] = logL'+(A) + C+, (34) 
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where are integration constants. (For sign-fixing: every argument of 
log(-) lias to be positive for large A > 0.) Now, WKB asymptotics like (31) 
also hold for A — > +cxd at fixed g, and we have painstakingly normalised 
earlier to ensure canonical large-A behaviours on all sides; therefore, any 
extra additive constants like are banned by (16), and (32) results. □ 

3.2 Explicit Wronskian identity 

A second solution of (1) besides ■'tpx, now recessive for q +e~"''/^oo, can be 
specified through the first conjugate equation, (6) ior i — 1, [15, §2.7] as 

*A(g) = V'SUe'^/^g). (35) 

Because the asymptotic form (19) for ipx and its analogues for 'ifx{q) and 
the g-derivatives all hold over the positive direction, the constant Wronskian 

of "ifx a-nd ipx, 1^(*a,'0a) ^aV'a ~ ^'x^Xi can be deduced explicitly, just 
from their q — ?> +cxo asymptotics. [15, Thm 21.1 (iii)] Under the canonical 
normalisations and recalling (22) for /3_i('y), it is [20, including Corrigendum] 

iy(*A,V'A) = 2ie''^/^e''^''-i(^^)/2 (A^7^2). (36) 

This basic Wronskian identity is the other essential piece of the frame- 
work: if we now use the exact- WKB identities (32) within the alternative 
form W^(*A, V'a) = *a(0)V'a(0) - *a(0)V'a(0), this then translates (36) into a 
functional relation between spectral determinants, which ior N ^2 reads 

e+'^/^DW+(e-'^A)Z^-(A) -e-^^/^D+(A)Z}W-(e-^^A) = 2ie''^^-^^^'^ . (37) 

3.3 Exact quantisation conditions 

At first glance, the original Wronskian identity (36), or its sibling (37), look 
grossly underdetermined: one equation for two unknown functions. This view 
however proves a delusion here, thanks to the structure formulae (14). Their 
first effect is to reduce quests for determinants to quests for their zeros (the 
eigenvalues): in other words, we are to seek exact quantisation conditions. 
The latter problem then becomes central to the whole framework; we now 
sketch our solution for it (as valid for > 2). [21] 
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To analyse the eigenvalues {Ek} of the Schrodinger operator H on the 
half-line, we may divide the functional relation (37) at A = —E^ by its first 
conjugate partner, to get 

2argi^W±(-e-'^£;)-(^^_i(i;) = ^[^+|±ifc|)] for k ^ . (38) 

The formulae (38) have the outer form of Bohr-Sommerfeld quantisation 
conditions (smooth left-hand-side functions of E, sampled at integers). Now 
they are exact, albeit not closed (the left-hand sides invoke the determinants 
of the first conjugate spectra, equally unknown). 

However, in either (Neumann / Dirichlet) sector — there is a complete 
decoupling between them — , the set of all quantisation conditions for the L 
conjugate problems plus the L structure formulae (now used to reconstruct 
each determinant from its spectrum), add up to a formally closed system: 



i-i 



logZ}[^+i]±(-e-i^Ejf^) -logD[^-i]±(-e+'^Ef ) - (-l)V/3-i(^?) (39) 



Ar-2 



2 2(Ar+2)J ' 



log D^^^iX) ^ Jnnjj2 fog(4'' + A) + | log{E^S + A) 



(40) 



{a>0} 

for £ e Z/LZ and k,K The selfconsistency brought by (40) makes 

this a fixed-point problem, with parities decoupled, hence abstractly of the 
form A4^{S'^} = £'^, where E*^ ("compound spectra") denote the disjoint 
unions of all the L conjugate spectra 8^^^^ of a given parity (Fig. 1), and M.^ 
are some spectra-to-spectra mappings. Note: as a completely general fact, 
concrete such forms and mappings are highly non- unique; we just hope to 
find some exphcit such form(s) of (39)-(40) that behave nicely. 

The equations (40) furthermore require asymptotically correct arguments 
8"^, for the K — )■ +oo convergence to take place. 

In the simplest case, all homogeneous potentials V{q) = q^: [19, 7] the 
compound spectra amount to the original ones 8^ C M+, mappings 
exist in the real form: {E'f. > 0} i-^ {E'^ > 0} {k ''™^), given by the formulae 

IA^{E'^)= k + l±^ forfcTdd' ^here 

A^iE") E arg(^; - e"^^ E") for i , E" > 0, (TV ^ 2), ^^^^ 
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Figure 1: (from [20]) The (odd) "compound spectrum" for quartic potentials 
V{q) = + V2q^, at V2 = —10 (left), (center), +10 (right). (In each square the 
L = 3 sets f M- are seen finitely cropped, labelled by their respective ^-values 0, 
1, 2, and rotated to fit their relative positions in the system (39).) 



and these have now been proven globally contracting. [2] Next, for mod- 
erately inhomogeneous potentials: we have only tested quartic and sextic 
cases numerically, [20, §3] and still - but now without proof to date - ob- 
served contracting maps. As long as this holds, the unknowns in (39)-(40) 
(the eigenvalues and the determinants) get uniquely determined by the fixed- 
point equations, and effectively reached as limits under iterations of the con- 
tracting maps; thus the fixed-point equations qualify as (exact) quantisation 
conditions. Only in some larger--u regions do our current schemes become nu- 
merically unstable - a non-elucidated issue. (Note: all numerical mappings 
have to be finite-dimensional as well, which entails some approximation.) 

Again for homogeneous potentials, moreover, (41) mysteriously coincide 
with the Bethe-Ansatz equations for certain solvable models of 2D statistical 
mechanics. This has fostered a new research area, the ODE/fM correspon- 
dence, relating our exact- WKB solvability with complete integrability in the 
modern sense. [7, 4, 16, 8] 

We may then briefly list ID Bohr-Sommerfeld quantisation conditions 
for the eigenvalues -E^, in the common general form A{Ek) = A; + |, as: 

- old exact-WKB setting, special: the Ek get specified exactly via an explicit 
function A{E), but only for exceptional systems (e.g., harmonic oscillator); 

- old general WKB approaches, semiclassical: the E^ get specified via an 
explicit function or expansion A{E), but only asymptotically for k — > -|-oo; 

- our general exact-WKB formalism: the E^ get specified exactly via a func- 
tion A{E) itself specified exactly, but implicitly through a fixed-point equa- 
tion which, while infinite-dimensional, is explicit. 
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3.4 Solving for the unknown functions tjj 




Figure 2: (from [22]) Exact-WKB calculations for an eigenfunction ip{q) (the 
ground state in a double- well potential). Dashed lines (with right-hand scale): 
for reference, the potential V{q) = q^—bq^ and its lowest eigenvalue Eq ^ —3.41014 
(input); vertical lines mark the two inner turning points. Continuous line (with 
left-hand scale): the eigenfunction ijj, canonically normalized as (31), drawn by 
numerical integration of (1) (NAG routine D02KEF). + : exact- WKB values ip{q) 
according to (32): meaning, at every q, a numerical iteration limit of the fixed- 
point system (39)-(40) for the determinant value D~{—Eq). The plot validates the 
exact- WKB method in extreme quantal regimes: a ground state, turning points 
(on the other hand, an unexplained numerical instability sets in beyond \q\ ~ 1.7). 



We focus on finding the canonical recessive solution tjjx for arbitrary A 
(once a particular solution is known, others trivially follow). 

Equations (32) have specified ipx^q) (and ip'xiq)) as particular spectral 
determinants over the half-line [q, +oo); now the latter are of the type fully 
resolved by the exact quantisation conditions above. Consequently, the un- 
known wavefunctions ip are computable as well in the exact approach (as 
numerical tests confirm: Fig. 2). [20, 22] 
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3.5 Toward ID quantum perturbation theory 

On potentials like V{q) = q^+gq"^, it is possible to monitor the 5^ — >■ 0+ limit of 
the exact- WKB framework, in spite of its extremely singular character (thus 
the degree A^, which is the most critical parameter in the framework, jumps 
from 4 to 2 at (? = 0). Results are largely governed by the singular behaviours 
of improper action integrals, like {q'^ + g^'^^^q'^ + Xy^"^ dq {g — )■ 0+) for 
the case V{q) = q^ + gq^. [23, 25] 



Appendix: derivation of (33) as trace formulae 

The identities (33) awaiting proof indeed take the form of trace formulae, 

MO) /MO) = Tr{H + X)-_' , MO)/MO) = Tr{H + X)'' , (42) 

where ip\ is the solution (31) of {H + X)ip = recessive for q +oo, and 
{H + A)^^ are the resolvent operators (Green's functions) on the half-hne q e 
[0, +00) with the corresponding boundary conditions at g = (— = Dirichlet, 
+ = Neumann). The resolvents are trace-class under our assumption N > 2. 

Note: the ■(/'-iiormalisations, so essential in the main text, inversely be- 
come irrelevant here; already in (42), any recessive solution cijjx will do. 

If now i>x{q) is the D*j™hiet solution of {H + A)^ = 0, also unique up to a 
factor, then the integral kernels G'^{q, q) of the Green's functions are known, 
in this ID setting, to be 

G^iq^q) = W^(^A,^J)-VJ(min(g,g))^,(max(g,g)) (g, q e M+), (43) 

where W{i/j, 0) '== ipcp' — (pip' denotes the Wronskian of {ip, (f)) (a constant 
when ('?/', 0) is a solution pair). As intermediate result, then, 

poo poo 

Tr{H + X)l'= G^{q,q)dq^W{M^^)-' i^^{q)M<l)M- (44) 
Jo Jo 

Next, the subtraction of tpx [{H + X)^)^] = from V'a [{H + A)'0a]' = 
(two obvious equalities) yields 

V'^V'A = [W-(V'^V'A)r; (45) 
this reveals an explicit primitive for the product V'a V'a, reducing (44) to 

Tr {H + X)l' = W{M ^t)~' ^(V^^ V^a) . (46) 
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But as g — )■ +00, using (17) and its derivatives, W{t/j^, — 
which tends to under N > 2, hence 

IV + X)l' = -W{^x, V'J)"' W{^1 ^x),=o. (47) 

Finally, this collapses to the identity pair (42) when W{ijjx,'i/j^) (constant) is 
itself expressed at g = and the boundary conditions for resp. ip\j 
applied. □ 
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